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Seiberg and Witten's proposed solution of TV = 2 SQCD with N c = 2 
and Np = 4 is known to conflict with instanton calculations in three dis- 
tinct ways. Here we show how to resolve all three discrepancies, simply by 
reparametrizing the elliptic curve in terms of quantities and u rather 
than r and u = (TrA 2 ). SL(2,Z) invariance of the curve is preserved. 
However, there is now an infinite ambiguity in the relation between and 
r and between u and ix, corresponding to an infinite number of unknown 



coefficients in the instanton expansion. Thus the reinterpreted curve (un- 
like the cases Np < 4) no longer determines the quantum modulus u as a 
function of the classical VEV a. 
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1. Review of discrepancies between the Seiberg-Witten solutions and instanton 
calculations 

The solutions of N = 2 supersymmetric QCD with gauge group SU(2) proposed by 
Seiberg and Witten [1] have recently been tested against first-principles instanton calcula- 
tions, up to the two-instanton level [2-4]. As reviewed below, these tests have resulted in 
perfect agreement when ATp, the number of quark hypermultiplet flavors, is < 2. On the 
other hand, interesting discrepancies have emerged for the cases N F = 3 [4] and N F = 4 
[3]. In this paper, we present our conjecture for how to resolve these discrepancies for 
N F = 4. Just like the conjectured resolution for N F = 3 [5], we will not actually need to 
modify the elliptic curve that governs the solution (Eq. (16.35) of [1]). Instead, we will 
change the interpretation of the parameters r and u that appear in this expression; they 
will have a different (and highly underdetermined) relation to physical observables of the 
microscopic theory than that proposed by Seiberg and Witten. In particular, in sharp 
contrast to the cases N F < 4, the reinterpreted curve for N F = 4 no longer determines the 
quantum modulus 

« = <Tr4 2 > , (1) 

as a function of the classical VEV a of the adjoint Higgs A = A a r a /2. 

The solutions given in [1] make precise predictions for all multi-instanton contribu- 
tions to the low-energy physics. The functional form of these contributions is tightly 
constrained by holomorphy, the U(1)r anomaly, renormalization group (RG) invariance, 
and the Matone relation [6-8] between the prepotential J 7 (a) and the modulus u(a). For 
N F > 0, additional constraints come from flavor symmetries, the decoupling limit for 
heavy flavors, and a discrete Z 2 symmetry forbidding odd-instanton contributions when 
any quark mass vanishes. As we showed in Refs. [2,3,9,10], all these constraints are built 
into the instanton calculus as well. 1 Therefore, in testing the Seiberg-Witten solutions 
against explicit instanton calculations, one can "mod out" these constraints, and focus on 
a basic subset of nontrivial numerical predictions extracted from the elliptic curves at each 

1 The all-instanton-orders proof of the Matone relation given in [10] was ostensibly limited to 
the case N F = 0. However, the reader can easily verify that the proof goes through for arbitrary 
numbers of (massless or massive) hypermultiplets. The key point here is that Eqs. (20)-(21) of 
[10] can be directly extended to N F > 0, see Eq. (7.20) of [3]. The Matone relation can also be 
understood as a Ward identity for (broken) superconformal invariance [11]. 
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order in the instanton expansion. These numbers are then compared against the results of 
a finite-dimensional integration over the (multi-) instanton super-moduli. 



N F 



parameter 
fixing 



first nontrivial prediction first nontrivial prediction 
of curve (>1 massless quark) of curve (all m q > 0 ) 



0 1-loop 

1 1-loop 

2 1-loop 



1-inst 



N/A 



2-inst 



2-inst 



2-inst 



2-inst 



3 l-loop+2-inst 

[Eq. (4)] 

4 1-loop 

+ all-even-insts 



none if all m q = 0; 
otherwise 4-inst 



4-inst 



3-inst 



3-inst 



[Eqs. (6), (42)] 



Table 1. For each N F < 4, the second column gives the order in the semiclassical 
expansion at which the input parameters and u (or and u for N F = 
4) in the Seiberg-Witten curves are fixed in terms of physical quantities in the 
microscopic theory. The third and fourth columns (respectively, at least one 
massless quark, and all quarks massive) then give the order in the semiclassical 
expansion at which a nontrivial numerical prediction is first made by the curve, 
which can be compared against an instanton calculation ("nontrivial" is defined 
in the text). To date, no numerical tests beyond the 2-instanton level have been 
performed, so that the models with N F = 0, 1, 2 are on a firmer footing than the 
ones with N F = 3,4. 

In formulating these tests, and evaluating their outcome, it is crucial to distinguish 
between what is input (i.e., the fixing of conventions) versus what is output (i.e., a definite 
numerical prediction) . Since this distinction will be key to resolving the discrepancies for 
N F = 4, we should first review the models with N F < 3 ; see Table 1 for a summary. The 
input parameters for the elliptic curves with 0 < N F < 3 are the dynamically generated 
scale A^ and the modulus it, as well as the N F quark masses. (Notation: henceforth 
we will denote by u the parameter in the Seiberg-Witten curves, and reserve the symbol 
u always to mean (Tr A 2 ).) The curve defines an implicit choice of regularization scheme 
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for AjJF. For any physical quantity the n-instanton contributions are proportional to 
(A^) n ( 4-iVF ) ; in order to test the specific numerical predictions of the curves against 
instantons, one needs to know how to relate to (say) the analogous quantity A^ in 
the Pauli-Villars (PV) scheme which is most natural for instanton calculations [12]. As 
explained by Finnell and Pouliot [13] and reviewed in Sec. 2 below, this "dictionary" is 
fixed at the 1-loop level, and reads: 

(A sw = 4(A^) 4 -^ , N F = 0, 1, 2, 3 . (2) 

Similarly, for 0 < N F < 2, a discrete symmetry on the u plane combined with a semiclas- 
sical analysis of the singularities gives simply [1] 

u = u , N F = 0,1,2 . (3) 

Equations (2)-(3) are the input; the predictions for J 7 (a) and u(a) first lie at the 1-instanton 
level for N F = 0, and at the 2-instanton level for N F = 1,2 due to the Z 2 symmetry 
mentioned above. 2 The tests for N F = 0 have been performed at the 1-instanton [13] 
and at the 2-instanton [9] level, and the tests for N F = 1,2 have been carried out at the 
2-instanton level [2-4], all with exact agreement. 

Next we consider the case N F = 3. Here Seiberg and Witten (in Sec. 14.1 of [1]) 
continue to assume the equality (3), even though it is no longer required by a symmetry 
argument. Indeed, for N F = 3, the Z^_^ F symmetry that acts on the complex u plane 
is trivial [1]; consequently Eq. (3) can be generalized to u = u — uq • (A| w ) 2 where uq is a 
numerical constant. (This possibility is noted by Seiberg and Witten, but not exploited.) 
A nonzero value for uq has been extracted from a 2-instanton calculation by Aoyama et 
al; they find [4] 

u = u-u 0 '(AD 2 , no = , N F = 3 (4) 

contradicting the naive assumption (3). It is natural to hypothesize [5] that the solution 
to the N F = 3 model is still given by the Seiberg- Witten curve (Eq. (14.5) of [1]), with 
the substitution (4) rather than (3). In sum, for N F = 3, both the 1-loop level and the 

2 When quark masses are nonzero the models with N F = 1,2 have 1-instanton contributions 
too; however, due to built-in heavy-quark decoupling properties, these contributions are pegged 
to the 1-instanton term for N F = 0 and so are not independent tests [2,3]. 
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2-instanton level should be considered as input due to the extra degree of freedom uq] the 
first testable predictions extracted from the curve then lie at the 3-instanton level (or, if 
one of the quarks is massless, at the 4-instanton level). 

In this paper we will focus on the particularly interesting case N F = 4. Here, unlike 
the previous cases, the proposed solution and the instanton calculation disagree at the level 
of the low-energy effective U(l) Lagrangian £ e ff [3]. In fact, this disagreement can already 
be seen in the massless model. In this case both the /3-function and the U(1)r anomaly 
vanish identically, so that the microscopic coupling #4 and ^-parameter 64 (which can no 
longer be rotated away) assemble into a single scale-independent holomorphic coupling 
r = + ^ . Furthermore £^ is simply the Lagrangian of a classical N = 2 free field 
theory; its overall normalization, r e ( ^, enters the BPS formulae for the dyon masses [1]. 
(The superscript (0) will denote the massless case.) In Section 10 of [1] Seiberg and Witten 
make the strong additional assumption that, thanks to the absence of a running coupling 
constant, the effective U(l) coupling equals the classical SU(2) coupling, 

*3 =r, (5) 

which implies no quantum (perturbative or non-perturbative) corrections to r. Instead, as 
we demonstrated in Ref. [3], a first- principles instanton calculation gives 3 

t3 = \F (0) "{a) = r + - J2 c ^ n . Q = ex P(^) • ( 6 ) 

77 n=0,2,4... 

In particular a nonzero 2-instanton contribution C2 = — 2^35 was calculated in Ref. [3]. 
The 1-loop perturbative constant cq, while not considered in Ref. [3], will turn out to be 
crucial to our proposed resolution. We calculate below that Co = 4 log 2 in the PV scheme. 4 
A second disagreement in the N F = 4 model arises when at least one of the four 
hypermultiplets has a nonzero mass, say rru ^ 0. In the double scaling limit defined by 
rri4 — > oc together with g 4 — > 0 in a specific way reviewed below, the heavy flavor decouples, 
and the model is supposed to flow to the N F = 3 theory. This requires the identification 

3 Note that we use q for the 1-instanton factor rather than for the 2-instanton factor as in [1]. 

4 If the classical exactness assumption (5) were correct, there would be no need to introduce 
a regularization scheme, since scheme dependence is a one-loop effect. As discussed in Sec. 3, 
although the iV> = 4 model is a finite theory, the microscopic SU(2) coupling g± must still be 
defined in a particular scheme; that is because the finiteness is due to cancellations between 
individually UV divergent graphs. 
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of the dynamical scale for N F = 3, with the parameters of the N F = 4 theory. As 
explained below, working in the PV scheme and using the dictionary (2), one obtains: 

A| w = 4Af v = 4m 4 exp(-87r 2 /(^ v ) 2 ) , (7) 

where g\ x is the microscopic PV coupling in the 4-flavor model. In contrast, the relation 
given in [1], 

A§ w = 64m 4 exp(-87r 2 /(^ w ) 2 ) , (8) 

involves a proportionality constant of 64 rather than 4. Again, if the quantum corrections 
to Eq. (5) were absent, there could be no distinction between the PV coupling constant 
g\ v and the "classical" #f w of Seiberg and Witten, hence no accounting for the factor of 
16 mismatch between Eqs. (7) and (8). 

Thirdly, even with the identification (8), the specific N F = 4 solution proposed in [1] 
flows to the uncorrected version of the N F = 3 model which fails to incorporate the shifted 
definition of u given in Eq. (4) . 

Our principal aim is to explain how to resolve all three of these N F = 4 discrepancies 
in a simple way, through a reinterpret at ion of the quantities r and u that enter into the 
massive Seiberg- Wit ten curve. In particular, rather than being modular forms of the 
microscopic SU(2) parameter r, the coefficients of the massive curve will be functions 
of the effective massless Z7(l) coupling defined by the all-even-instanton series (6). 
(Obviously the factor of 16 between Eqs. (7) and (8) will be automatically accounted for 
by the exponentiation of c 0 = 4 log 2.) This redefinition preserves the important £X(2, Z) 
invariance of the elliptic curve, as well as the stringent residue condition described in Sec. 17 
of [1]. However, it also introduces an infinite ambiguity into the solution, parametrized by 
the infinite number of as-yet-undetermined numerical coefficients C4, c 6 , . . . , in Eq. (6) (a 
similar series, Eq. (42), relates u and u). These issues are discussed in Sec. 3 below, which 
is devoted to the case N F = 4, and also contains our conclusions. But first, in Sec. 2, we 
lay some necessary groundwork in the cases N F < 4. In particular we review Weinberg's 
matching prescription between high- and low-energy gauge theories [14], as this formalism 
lies at the heart of the physics. In the process we will also specify the perturbative, one- 
loop contributions to the prepotential for N F < 4, correcting some incomplete expressions 
in the literature. 

The main result of this paper is the reinterpreted N F = 4 massive curve which agrees 
with all available perturbative and multi-instanton calculations. Although the input pa- 
rameters of this curve receive contributions from all even numbers of instantons, the curve 
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does contain definite numerical predictions at the 3-instanton level (assuming nonvanish- 
ing quark masses) that can, in principle, be tested against a semiclassical calculation. We 
anticipate that similar reinterpretations need to be made in the general class of N = 2 
models with gauge group SU(N C ) and Np = 2iV c for which the /3-function vanishes. 

2. Weinberg's matching prescription and the Seiberg-Witten regularization 
scheme for Np < 4 

The physics of N = 2 SQCD utilizes — in two distinct but equally important ways — a 
matching prescription between a "high-energy" and a "low-energy" gauge theory. On the 
one hand, the 577(2) gauge group spontaneously breaks down to U(l) as the adjoint Higgs 
A acquires a complex VEV ( A) = ar 3 . For energy scales E < M w = y/S\a\, the dynamics 
is governed by a nonrenormalizable Wilsonian effective action with 17(1) gauge invariance, 
formally obtained by integrating out the heavy quanta. On the other hand, one also needs 
to understand the RG decoupling mentioned earlier, whereby the SU(2) theory with N F 
flavors of quark hypermultiplets flows to the SU(2) theory with N F - 1 flavors, in the limit 
that one of the quarks becomes infinitely massive. 

Both types of matching may be accomplished with the use of Weinberg's one-loop 
formula [14] 

Here /x is the characteristic momentum scale of the light fields; £ L e(aO and # H e(aO are 
renormalized gauge couplings of the low- and high-energy theories, respectively; and A(/z) 
is a finite correction coming from one-loop contributions of heavy particles to the gauge 
self-energy E. Thus g HE (/x) is extracted from the complete set of one-loop contributions to 
E, whereas for £ LE (/x) only light particles are permitted on the external and internal legs. 
For example, in the MS scheme, £ H e(aO and £ le (a0 are defined in D dimensions in terms 
of the bare couplings # H eb an d 9leb m the standard way: 

9he B ^ 2 " D/2 = 9he(^) - &he glM {^ZTl + 57e - 5 lo § 47r ) ' ( 10 ) 

g^ B v 2 ~ D/2 = 0 LE (/x)- h*glM {jy^l + ^ 7e ~ ^ log4?r ) 5 ^ 

where 6 HE and 6 LE are the one-loop coefficients of the corresponding /3-functions. Note 
that when the low-energy theory is supersymmetric pure 17(1) gauge theory, # LE receives 
no perturbative corrections and is scale-independent in the limit D — > 4. 
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The quantity A has the generic form 

X(fi) = C 1 + C 2 log^ + C3log^+C4log^ , (12) 

jjj [Jj jJj 

where m v , rrif and m s are the masses of the heavy vector, fermion and scalar particles 
that have been integrated out. The Q are group-theoretic constants that are tabulated 
by Hall (see Appendix 1 of [15], in which an error in [14] is corrected). In particular C± 
and C 2 appear only when the heavy vector particles have been integrated out; similarly, 
C 3 appears when there are heavy fermions, and C4 corresponds to the heavy scalars. 

The original calculations of [14-15] were performed in the dimensional regularization 
with MS scheme (DREG). However, they can be easily translated into the supersymmetry- 
preserving dimensional reduction with MS scheme (DRED) [16], or into the Pauli-Villars 
(PV) scheme. The one-loop relations between the coupling constants in different schemes 
for the SU(N C ) gauge group can be found for example in [17,13]: 

L_ = _J_ = _J_ + J_ N (13) 

independently of Np. It turns out that for the case of spontaneous symmetry breaking 
[SU(2),N F ] — ► [17(1), 0] considered in detail presently, 

In contrast, in the case of the heavy flavor decoupling [ST/ (2), N F ] — > [SU(2), N F - 1], 
vector particles do not decouple and Ci = 0 by definition. 

Now let us apply this formalism, in turn, to these two cases of interest. 

1. Spontaneous symmetry breaking 

We consider the Coulomb branch of N = 2 SQCD, in which only the adjoint Higgs 
A = A a r a /2 acquires a VEV, say in the r 3 direction: (A) = ar 3 . (We have adopted 
here the VEV normalization conventions of [1]; the translation formulae to the original 
normalization of [18] used in our previous work are assembled in Appendix A.) The SU{2) 
component of the adjoint N = 2 supermultiplet that is parallel to the VEV then remains 
massless, while the components oc r 1 or r 2 acquire a mass M w = V8\a\. In addition there 
are 2N F quark multiplets Qf and Q/, / = 1, ■ • ■ , JV>, in the fundamental representation 
of the gauge group. The '1' and '2' color components of these multiplets acquire masses 
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I \Fla + rrif \ and | \/2a — m/|, respectively, as can be seen from a tree-level examination of 
the N = 2 invariant superpotential [1] 

N F 

W = J2V2Q f ^Q f + m f Q f Qf . (15) 
/=i 

Here <& is the N = 1 adjoint chiral superfield whose lowest component is A; color indices 
are suppressed. 

In the PV or DRED schemes one then has 

am = cu log ^ + c„„ d e (log h^uml + log h^M) , (16) 

where C adj and C fund are numerical constants. Setting /x = for simplicity and extract- 
ing C fund = 1/167T 2 from Ref. [15], one finds 



1 N F 



2a 2 - m 2 



8a 2 



(17) 



so that Eq. (9) becomes 

1 i-ioop 4 — Np 1 /M ff \ JV f 1 ^ , 2 2| 



(18) 



In this expression # eff = # LE denotes the effective E/(l) coupling constant, A Nf is the PV 
or equivalently the DRED dynamical scale (we drop the PV superscript henceforth) 

A b Nr = /x 6 exp[-87r 2 /^ v(/x)] = ^exp[-87r 2 /^ RED (^)] , (19) 

and we have used the fact that the (negated) coefficient of the /3-function for these models 
is b = 2N C - N F = 4- N F . We make the following comments: 

(i) Equations (18)-(19) extend to N F > 0 the case of N F = 0 considered by Finnell 
and Pouliot, who obtained simply [13] 



1 l-loop 1 = J_ ] (M w \ 

<4 9lv{M w ) ^ ° g V A 0 ) 

This is referred to as the absence of threshold corrections. 



log == . (20) 
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(ii) As usual in a supersymmetric chiral theory [19], l/g^g * s the imaginary part of a 
holomorphic complexified coupling constant 

TefT = "2- + — • 21 

This implies that Eq. (18) may be analytically continued away from the imaginary axis, 
as follows: 

r eff (iM(4-N F )log(^j + ^log2 - ]og(l-mJ/2a a ) . (22) 

£ Heavy flavor decoupling 

Next we consider the case in which a single hypermultiplet flavor becomes very heavy 
(say, Tftjsfp — » oo) and decouples from the spectrum, leaving behind N = 2 SQCD with one 
fewer flavor. Choosing fi = m NF we find that \{m NF ) = 0 as follows from Eq. (12) with 
Ci = C 2 = 0; consequently 

9N 2 F (m NF ) = 9xl-i{™>N F ) (23) 

in the PV or DRED schemes. Rewriting this relation in an RG-invariant way using (19), 
one obtains 

m^.AjT/' = A^r 1} . (24) 

The appropriate double scaling limit is therefore defined by mjsr F — >• oc and Aiv F — > 0 with 
the product on the left-hand side of (24) being held fixed. 

In the remainder of this section we will relate the dynamical scales AJJF that appear 
in the Seiberg-Witten elliptic curves, on the one hand, to the PV or DRED scales (19), 
on the other hand. This is done by comparing Eqs. (22)-(24) to the explicit solutions 
proposed in [1]. 

3. Relating the Seiberg-Witten scheme to the PV or DRED schemes 
The Seiberg-Witten elliptic curve for the N F = 0 theory is [1] 

y 2 = x 2 (x-u) + J(Ag w ) 4 s. (25) 

The curve defines a dynamical scale A^ w in a particular scheme, the "Seiberg-Witten 
scheme" ; the superscript SW is introduced to distinguish it from A 0 in the P V or DRED 



10 



schemes. The correspondence between the SW and the PV schemes for Np = 0 has 
been examined by Finnell and Pouliot [13], using Weinberg's matching formula. From the 
elliptic curve (25), they calculate r e ff ~ (2i/ir) log (8u/(Aq w ) 2 ) valid in the semiclassical 
regime, u ~ 2a 2 . Comparing this result with Eq. (22) with Np = 0 then yields [13] 

(AH 4 = 4 A£ • (26) 

While this derivation was ostensibly performed at the one-loop level, it is well known that 
such relations between A's defined in different renormalization schemes are actually one- 
loop exact [20], regardless of the presence of supersymmetry. (In contrast, the definition 
(19) of the dynamical scale itself is one-loop exact under the RG of the Wilsonian effective 
action, but only because of supersymmetry [19,21].) 

Next we consider the cases 0 < N F < 4. Let us introduce the symmetric polynomials 
in the masses: 

N F N F N F 

= 1 , = J>? , = J2 m\m) , • • • , M<g'> = J] ™3 ■ (27) 

*=1 i<j j=l 

The curves are then given by [l] 5 

f iNF) = x 2 (x-«) + ^0^(A^) 4 -^x-^(A-) 8 - 2 ^ J2 MF'Hx-ii)*'- 1 -' 

6=0 

(28) 

Seiberg and Witten simply equate u = u for Np = 1,2,3 but, as reviewed earlier, for 
Np = 3 this should be corrected to Eq. (4). It is easily checked that in the decoupling 
limit Tfi]s[ F — Y oc, one obtains the desired result yf Nr) —> yf Nr _ 1) if and only if one makes 
the identification 

m NF ■ (AZ) 4 ~ Nf = (A^_i) 4 - (iW) • (29) 



5 A technical aside: In the absence of quark masses we have identified alternative curves for 
both Nf = 2 and Nf = 3 with the desired singularity structure. They are y 2 — x 2 (x — u) — 
^(A2 W ) 4 (ar — u/9) and y 2 = x 2 {x — u) — ^(A| w ) 2 u 2 , respectively. Consistent with the physical 
arguments in [1], in the former case the discriminant A (it) correctly factors into a product of two 
double roots, while in the latter case A(u) is the product of a simple root with a quartic root. 
However, in each case A(u) has the wrong singularity structure once mass terms are added, which 
can be seen most easily in the case that all masses are set equal. 
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Notice that this is the same recursion relation as for the PV and DRED schemes, Eq. (24). 
From the Np = 0 "boundary condition" (26), one immediately derives the dictionary 
between schemes given in Eq. (2) above. 

Before passing to the case Np = 4, we note that the perturbative, one-loop structure 
of the effective 17(1) coupling contained in Eq. (22) has rarely appeared correctly in the 
literature. This expression has recently been confirmed by [22]; rather than invoke the 
Weinberg prescription, these authors extract the result directly from the elliptic curves, 
for arbitrary N c . 

3. Resolving the discrepancies in the Np = 4 solution 

We now turn to the interesting case N F = 4. In this model the (3- function vanishes 
and no dynamical scale is generated. Note that it is trivial to extend the RG matching 
relation (24) to this case. The relation (23) still holds when ji = m 4 — » oc; consequently 



using Eq. (19). Here, and henceforth, g^) = 9a is the scale-independent microscopic cou- 
pling in the PV or DRED scheme 6 ; it combines with 9 4 to form a single scale-independent 
PV or DRED holomorphic parameter 



Before turning to the elliptic curve, let us revisit the three discrepancies with instanton 
physics that we wish to address. The first is the discrepancy between Eqs. (5) and (6) 
which relate and r. (As above, we will use the superscript (o) to denote the massless 
case.) The one-loop constant Co in Eq. (6) was not considered in Ref. [3]. In light of the 

6 One may ask why, in this finite theory, it is nevertheless necessary to specify a scheme. As 
mentioned earlier, this is because the finiteness is due to an "oc minus oc" cancellation between 
divergent diagrams, which is intrinsically ill defined. For this reason we cannot simply construct 
physical quantities directly from the bare tree- level coupling gs, as would be natural to do if all 
individual graphs converged. That gB cannot be a scheme-independent physical quantity can be 
seen from Eq. (10), which applies equally to the DRED and DREG schemes, and which would 
then imply for N F = 4: # D red(aO = #dreg(m) = 9b /J> 2 ~ D/2 in contradiction to Eq. (13). 



m 4 exp(-87r 2 /#!) = A 3 



(30) 



8iri 9 4 



(31) 
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preceding discussion, we can now read off the value cq = 4 log 2 from Eq. (22), which is 
easily extended to the case N F = 4: 

4' * 4 

r eff = \T"{a) = r + ^log2 - JL £ log ( 1 - m)/2a 2 ) + (9(g) . (32) 

An independent derivation of this 1-loop relation (which tests our proposed reinterpret at ion 
of the massive curve) is discussed below. 

The second and third discrepancies involve properties of the massive N F = 4 curve 
(Sec. 16.3 of [1]). Consider the RG decoupling property when the quark mass tua grows 
large. In the double scaling limit the curve indeed collapses to the N F = 3 curve (28); 
however, as reviewed in Appendix B, this limiting behavior requires the identification [1] 

64m 4 exp(-87r 2 /(^ w ) 2 ) = A| w = 4A 3 (33) 

where the second equality follows from Eq. (2). This apparently contradicts Eq. (30) above; 
more precisely it means that the coupling #f w used in [1] cannot be the expected DRED 
quantity g±. And thirdly, the N F = 3 curve that one flows to in this way has u = u rather 
than the shifted definition (4). 

As the reader can anticipate, we will posit that the parameter r that appears perva- 
sively in the Seiberg-Witten curve for N F = 4 should really be identified with the effective 
massless Z7(l) coupling r^, Eq. (6), rather than with the microscopic SU(2) coupling r, 
Eq. (31). By definition, this reinterpretation resolves the first of the three discrepancies. 
Pleasingly it also resolves the second discrepancy: the constant factor c 0 = 4 log 2, when 
exponentiated, precisely compensates for the factor of 16 mismatch between Eqs. (33) and 
(30). (A more stringent test of our proposal is discussed below.) Finally, the third discrep- 
ancy will be resolved by altering the relation proposed in [1] between u and u\ however, 
there is an infinite ambiguity in this procedure that we do not know how to eliminate. 

We now review the Seiberg-Witten curve (Sec. 16 of [1]), starting with the massless 

case: 

{y (0) f = x 3 - \g 2 (r syf ) xu 2 - \g 3 (r syf ) u 3 = W^W^W™ , Wf > =x- ei (r sw ) u . 

(34) 

Here g 2 and <7 3 are rescaled Eisenstein series. The cubic roots may be denned in terms 
of ^-functions; they have the semiclassical expansion 

ei(r) = | + 16g 2 + 16g 4 + 0{q G ) , e 2 (r) = -± - 8q - 8q 2 - 32q 3 - 8q 4 + 0(q 5 ) , (35) 
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with e 3 = — ei — e 2 . As noted in [1], the ei are not strictly speaking modular forms of 
SL(2,Z). Rather, they are weight-two modular forms of three different conjugate sub- 
groups of SX(2, 2Z)\ under the action of the full group they permute amongst themselves. 
The curve (34) is well known in the math literature [23]. It is designed so that if the VEVs 
a and dp are extracted in the standard way as periods of the curve [1], 

da y/2 f dx dan a/2 f dx ( . 

then one has simply 

a = v^/2 5 ^ = r sw a . (37) 

As expected, these are the defining equations of a classical free field theory, with !F i0) (a) = 
t sw a 2 and u = 2a 2 . 

Seiberg and Witt en make the two further assumptions r sw = r and u = u. The first 
of these assumptions contradicts Eq. (6); instead, we will assume r sw = as explained 
above. As pointed out in [5], this, in turn, invalidates the second assumption as well; 
instead, one must take 

u = «■ (dr^/dr)' 1 . (38) 

This latter redefinition is specific to the massless model; it follows directly from the in- 
stanton version [24,10,5] of Matone's relation, which for four massless flavors reads, 

u — 2mq— — , (39) 

combined with Eq. (37) . 

We turn finally to the massive curve. Setting eij — e% — ej, one has [1]: 

y 2 = W x W 2 W 3 + e 12 e 23 e 31 (Wtf&s + W 2 T 2 e 31 + W 3 T 3 e 12 ) - e 2 12 e 2 23 e 2 31 N . (40) 

Here Wi = W- 0) - e 2 i?, where i?, N and are symmetric polynomials in the four masses: 

R = 1^2™* ' N = Ye ^ m ? m i m fc-^S m M + ^S m ^ 

» i>j>k i^j i , s 

i>j i i i>j i 

with T 3 = -Ti - T 2 . Seiberg and Witten observe that the coefficients of this curve are 
£X(2, Z) modular invariants, provided that the aforementioned permutations amongst the 
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ei are accompanied by the same permutations acting on the (this is referred to as SO (8) 
triality). They view this as strong circumstantial evidence that the dyon spectrum itself 
is SL(2,Z) invariant. 

In restoring agreement with instanton calculations we need not tamper with these 
important SL(2,Z) properties. Rather, we will reinterpret the intrinsic parameters of 
the curve. In particular we will continue to take = e^(r e ( ^) where r e ( g? is given by the 
instanton series (6). It remains only to relate the parameter u (which enters through the 
Wi) to the physical quantum modulus u = (TV A 2 ). Dimensional analysis, 0(4) symmetry, 
smoothness in the masses, and the above-given massless limit suggest the following generic 
relation: 

t/r (0) _i 

u = u- + R- J2 a ™<l n ' S = exp(i7rr) . (42) 

T n=0,2,4-. 

(This expression is to be compared with the Seiberg-Witten proposal [1] 

u = u - \ei{r)R (43) 

which is already faulty in the massless limit.) The absence of odd instanton contributions 
is due to the discrete Z 2 symmetry discussed earlier, the mass parameter R being even 
under this symmetry. The numerical coefficients a n in (42) may be constrained using a 
variety of physics considerations explained in Appendix B. We find that a$ = —1/3, and 
a 2 = 37/(3 3 2 5 ) ; the latter value disagrees again with the proposed relation (43), but is 
necessary to recapture the shifted definition (4) of u in the decoupling limit. 
We conclude with the following comments: 

(i) The higher-instanton contributions a n with n = 4, 6, • • • remain completely un- 
determined, as do the constants c n with n = 4, 6, • • • in the relation (6) between r e ( ^ 
and r. When masses are incorporated into Eq. (6), Matone's relation [6,7,10] may give 
an interesting correspondence between the two series (as it already does in the massless 
case [5]). 

(ii) Our lack of complete knowledge of the relation between u and u does not actu- 
ally affect the low-energy effective Z7(l) Lagrangian. This is because a and an are still 
determined by Eqs. (36) (with y instead of y (0) in the massive case). Both sides of these 
equalities involve the (independent) variable u rather than the unknown (dependent) vari- 
able it. The former is then eliminated in favor of a, giving a prepotential in which neither 
u nor u appears: T = T{a, w(o); = F(a;{mi}). In other words, so long as the 
classical VEV a is considered the independent variable (and not the quantum modulus 
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u = (Tr<fr 2 ) on which a depends in a presently unknown way, and vice versa), the prepo- 
tential is "known" as a function of (albeit not as a function of the microscopic r). In 
contrast, for Np < 4, both a and T are known functions of u as well. 

(iii) As a stringent test of our proposed redefinitions ei = ei(r^), we have in fact 
constructed F{a) as outlined in (ii), and verified that the right-hand side of Eq. (32) is 
indeed reproduced by the curve (paralleling Ref. [22]). 

(iv) Finally we reiterate the point that the massive Np = 4 curve is still SL(2,Z) 
invariant, but only in terms of the rather non-intuitive quantity r e ( g? (the effective coupling 
in the massless theory) rather than the microscopic coupling r. While the relation (6) 
between these parameters is currently unknown beyond the 2-instanton level, it would be 
pleasing if, in the end, the model turned out to be modular invariant in terms of r as well. 
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Appendix A. Note on conventions 

Here we provide the dictionary between the normalization conventions of Ref. [1] which 
we have adopted in this paper, and those of Ref. [18] which we used in previous work [9,3]. 
The original unrescaled VEV definitions of [18] are: 



(A) = vr 3 /2, (A D ) = v D r 3 /2, v D = 

87Ti 0eff 0 8V D 0 9 2 J-(V) 

T~eS = ~2~ H = — = ^ — » 

gl s 7T dv dv 2 

The VEV normalizations of [1] adopted in this paper are: 

a = 2V , a D = y d = 2~fa- 

dan t d 2 F(a) 
Teff " da ~ 5 da 2 

v = V2u+--- , a = |V2w+--- 
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dv 



The prepotential for Np < 4 massless flavors in the two different normalizations reads: 

^MR)^' 1 - £ F n NF) ( ^ ) « 2 

77 „=2,4,6 ) ... V G / 

(a \ n(4—N F ) 

Using Eq. (2), the relation between the instanton coefficients is determined to be: 

jr(JV» _ 2 n ( 6 - A? >)- 2 ^(^iO 

Appendix B. Constraints on the relation between u and u 

In this appendix we constrain some of the a priori unknown numerical constants a n 
that relate u to u as per Eq. (42). We use the following two considerations: 

(i) Closely following Sec. 16.3 of [1], we first consider the illuminating special case 
(mi,m 2 . ra 3 . ra 4 ) = (0,0, m,m). Then JV = = 0 and one has simply y 2 = WiW 2 W 3 
with Wi = x — eiii — e 2 m 2 . As always, the critical points on the quantum moduli space are 
the values of u for which two of the three x roots coincide. In the present instance this 
means e^u + e 2 m 2 = eju + e 2 m 2 , or equivalently since e i — 0 : 

u = {eim 2 ,e 2 m 2 ,e 3 m 2 } ~ {|m 2 , (-J - 8q)m 2 , (-J + 8g)m 2 } (B.l) 

using Eq. (35). We have dropped terms of order g 2 as these do not survive the double 
scaling limit. Now consider the decoupling limit m — » oc. On the one hand, one expects 
a perturbative singularity at u ~ 2a 2 ~ m 2 which corresponds, physically, to two quark 
multiplets becoming massless (see Eq. (32)). On the other hand, the leftover model after 
the two heavy flavors decouple is the massless N F = 2 theory, which has singularities at 
u = ±|(A1 W ) 2 (see Eq. (28)). In sum, 

u ~ {m 2 ,-|(Ai w ) 2 ,i(An 2 }. (B.2) 

Equating (B.2) with (B.l) forces a 0 = -1/3; and furthermore 8gg m 2 = J(A| W ) 2 which 
is precisely consistent with the DRED recursion relations (29), (30) and (33) (with #f w = 
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(ii) Returning to the case of four generic masses, we next consider the double scaling 
limit rri4 — > oo, — > 0 with the product q^ fixed at A| w /64 as per Eq. (33). Naively, 
we would like the curve (40) to collapse to the N F = 3 curve defined by Eqs. (28) and (4). 
Instead, in this limit the right-hand side of (40) diverges badly. However we can exploit 
the fact that the variable x is just a dummy of integration (see Eq. (36)); an appropriate 
shift in x eliminates this divergence and guarantees a smooth RG limit. Accordingly we 
let 

x -> x + (l + (±-l*2)q 2 + 0(q*))R - {\ + G{q*))u (B.3) 

which obscures the SL(2,Z) properties of the curve, but preserves the Z 2 properties. 
The explicit factors in (B.3) have the following genesis: the constant 2/9 eliminates the 
large-mass divergence; the constant 1/3 guarantees that the cubic terms in x and u will 
have precisely the form x 2 (x — u) dictated by Eq. (28); and finally the factor (1/32 — a 2 /3) 
eliminates the (A| w ) 3 mim2TO3 term that otherwise generically appears, again in order 
to harmonize with Eq. (28). (This term is not forbidden by any symmetry; its absence 
from the Np = 3 curve (28) is convention dependent, as is the form of the cubic term.) 
Substituting Eq. (B.3) into the Np = 4 curve (40) and taking the double scaling limit then 
indeed reproduces the N F = 3 curve (28), provided that 

-(l + 32a 2 )2- 10 = u 0 . (B.4) 

From the value of u 0 quoted in Eq. (4) one deduces a 2 = 37/(3 3 2 5 ) . 
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